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Abstract 

We reconsider the virial theorem in the presence of a positive cosmological 
constant A. Assuming steady state, we derive an inequality of the form p > 
A{K/AitGn) for the mean density p of the astrophysical object. The parameter 
A depends only on the shape of the object. With a minimum at A sp h ere = 
2, its value can increase by several orders of magnitude as the shape of the 
object deviates from a spherically symmetric one. This, among others, indicates 
that flattened matter distributions like e.g. clusters or superclusters, with low 
density, cannot be in gravitational equilibrium. 

Due to its wide range of application as well as its generality, the virial theorem 
plays an important role in astrophysics. To derive the virial theorem, one requires only 
the collisionless Boltzmann equation. Assuming steady state, one of the important 
applications of the virial theorem is to deduce the mean density of astrophysical 
objects like galaxies, clusters and superclusters by observing velocities of a 'test- 
body' around them. It is clear that for conglomeration of matter, spread over large 
enough scales, the Hubble expansion of the universe will, in principle, oppose the 
gravitational equilibrium. It is also known that a positive cosmological constant 
A > accelerates this expansion (for a review on the cosmological constant and 
its problems, see |L|). If, as recent measurements seem to indicate 0, a positive 
cosmological constant enters the Einstein's equations, the resulting space-time in the 
Newtonian limit (called Newton-Hooke space 0) will inherit the expansion due to 
the A-term in the form of an external force. In this limit of the Einstein's equations 
we can rederive the virial theorem and evaluate the conditions under which a steady 
state for a collection of matter is reached when, as it is the case here, we have two 
opposing forces: the attractive Newtonian force and the repulsive external A-force. In 
exploring the astrophysical significance of A > 0, we will make use of the revised virial 
theorem (other approaches to eventual astrophysical effects of A have been discussed 
in HQ and [§). 

Before going into medias res of the virial theorem, we mention first some salient 
features of the Newtonian limit itself, including now the cosmological constant H. 



1 



This limit is, in many respects, quite different from the A = case. The Newto- 
nian limit of general relativity is given in the form of a Poisson equation for the 
gravitational potential $ which has to satisfy 



|$| < 1 



In our case the Poisson equation reads 



A$ = (4irG N )p - A 



(2) 



with Gn being the Newton's constant. Trivially, to solve (|2|) we need some boundary 
conditions. In case of A = 0, the latter is taken to be of Dirichlet type and is put 
at an infinite distance R i.e. $|k=oo = const. This is consistent since any Newtonian 
potential will fall off at large distances as 1/r. For A / there is, besides the 
equation yet another source of information about the Newtonian limit. This is the 
Schwarzschild solution for of spherically symmetric bodies. Via goo — ~~ (1 + 2$), 
we can then infer from this solution the potential of a point-like and spherically 
symmetric body. It reads 

$(r) = -Ar 2 (3) 

Combining this with equation ([]]), we obtain two constraints of the validity region of 
the Newtonian limit, namely 
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Note that in case of A = neither M max nor R max exist. On account of R max > we 
cannot anymore put the Dirichlet boundary condition at infinity if we want to have 
a consistent Newtonian limit. We have to rather choose some finite distance R where 
we can set = 0. The solution of the Poisson equation (Q) with this boundary 
condition can then be written as 
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andr< = min(|cc|, \x'\), r> = max(|cc|, \x'\). Hence, there are now two effects of the A- 
term in the potential: a direct term in the form of an external potential — (1/ 6) A| | 2 
and an indirect one, involving G'(x,x'), which appears because we were forced to 
choose the boundary condition at a finite distance. If R is some fraction of R ma x 
(the only distinguished large radius available), but still quite sizeable, the indirect 
contribution to $ will be suppressed by powers R~ n and, in the first approximation 
we can neglect its appearance in the virial theorem. 

Having discussed the Newtonian limit to some extent, we derive next the 
virial theorem including the external — (1 /6)A|a;| 2 potential. We could start with the 
collisionless Boltzmann equation, but for simplicity we will follow here a more pedes- 
trian approach by considering a discrete collection of massive objects enumerated by 
a, j3 = 1, 2, 3, N. The inertia tensor is then given by 



N 

k 
a=l 



Differentiating this tensor twice with respect to time, using Newton's equation of 
motion for x° and defining the kinetic energy tensor Kjk and the corresponding 
potential energy tensor Wjk by 



1 A 

x ^ • a • ex. 
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Wjk = W kj = --G N 2 ™ a m p J_ (9) 



we readily arrive at 



d 2 Ijk . ~~ - 



ijf2 AK 3 k + 2Wjk + -M, k (10) 

which is the virial theorem. It is worthwhile stressing that given Kjk and Wjk the 
presence of A makes the virial theorem a differential equation for Ijk- However, 
lacking a priori this information, one can proceed as follows. First we take the trace 



of equation ([10]) and define K = Tr(K jk ), W = Tr(W jk ), I = Tr(I jk ). Next we 

dt 2 



assume a steady state i.e. = 0. Then equation (pi]) becomes 



2K + W + -AI = (11) 
3 

Had we started with the collisionless Boltzmann equation, the final result would have 
been the same. It is now convenient to pass to continuous systems for which 



I = dxp(x)\x\' 

K = - [ dxp(x)v 
2 Jv 



W = ~! dxp(x)$ N (x) (12) 
2 Jv 
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where v is the averaged velocity and <3>jv(£c) is the Newtonian part of the solution in 
(H) (proportional to Gjy neglecting the third term involving G'(x, x')). Since $at < 0, 
it follows that W < 0. On the other hand K is positive definite. This leads to the 
basic inequality of gravitational equilibrium with a non-zero A 

--AI+\W\>0 (13) 
o 

Note that this inequality is trivially satisfied if A = 0. Any astrophysical object 
which does not obey this inequality, cannot be in a gravitational equilibrium. Clearly, 
equation (|T3|) is a relation between the density of the body and its shape (geometry), 
the latter encoded in the integrals over the volume of the body. One way to exploit 
(|l3| ) would be to model the density for the object under consideration (as it is often 
done in astrophysics) by assuming a functional form p = p(x; Aj) where the Aj's 
are some parameters (often of dimension of mass or length). Then the inequality 
(p~3| ) would convert into a relation among these parameters. It is, however, more 
transparent to work with a constant mean density p. Defining then 

I = pl\ w = - l - P 2 G N W (14) 

we can recast the inequality (|13"D into the form 

P > Ap VSkC 

16tt I' A 
A = ~W>' Pvac = 8^ (15) 

To search for the effects of A, inequality ([15]) has an appealing form as it is solely 
expressed through the mean density and a parameter A which depends only on the ge- 
ometry. It is obvious that the larger the value of A is, the more will be the importance 
which we can attribute to A with regard to gravitational equilibrium. 

To investigate the possible sizes of A we consider the volume of the body 
bounded by smooth surfaces of second degree. With the triaxial ellipsoid many dif- 
ferent volumes can then be approximately modeled. For the sphere of radius Rq we 
obtain W = (32/15)n 2 R 5 and V = (4/5)^ resulting into 

^sphere = 2 (16) 

As we will see below, this seems to be the smallest possible value of A. For an 
ellipsoid with three axes a, b and c the integrals give W = (16/15)7r 2 a6cJa 2 and 
/' = (A/15)irabc(a 2 + b 2 + c 2 ) where the parameter J depends on the type of the 
ellipsoid i.e. on the remnant rotational symmetry. The parameter A for an ellipsoid 
has the general form 

4 a 2 _|_ jj 2 _|_ c 2 

^ellipsoid = ~ J - ^ (17) 



4 



We quote now the parameters J's for the three different types of ellipsoids 
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where F(<p, k) is the elliptic integral of the first kind. The corresponding three different 
parameters for A can be calculated to be 
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(19) 



To appreciate the final results (eq. ([l5|) together with (|19D ), let us take an example 
of a flattened prolate ellipsoid with (c/a) = 10(10 2 ). The resulting A pro i at e is then 
10 3 (10 7 ). Given first the fact that the nowadays preferred value of p vac is (0.7— 0.8)p cr i t 
(Pcrit = 1-789 x lCT 29 /iQgcm~ 3 , 0.6 < h < 0.8) and secondly that in comparison to p vac 
the density of a typical elliptic galaxy is (0.1 — 2) x 10~ 26 gcm -3 , we can safely state 
that inequality (|T5|), and hence the effect of A in considering gravitational equilibrium, 
is of astrophysical relevance. The type of astrophysical objects (galaxies, clusters and 
superclusters) we are investigating may or may not be in gravitational equilibrium. 
In both cases, the above results seem relevant. In the first case, with a fixed density 
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p, we would expect that the deviation of the rotational symmetry is indeed small. In 
trying to reach the equilibrium, the objects will reduce the ratio of their axes. This 
is a direct consequence of the inequality (15) and the values of A given in ([16]) and 



fll9f) . As such, this is also a consequence of A > 0. If, in contrast the inequality 
(|T5| ) is not satisfied, the gravitational steady state is not reached and we should be 
careful in extracting the densities from the virial theorem in such situations. First, 
for low density objects the effect of A should be included in the virial theorem (eq. 



11)). Furthermore, the inequality (|15|) can always serve as a check if the assumed 
equilibrium indeed holds. 

Let us now consider objects which are most likely not in a steady state. These 
are superclusters which appear mostly in a flattened form and have low densities. 
Even for richly populated superclusters with p ~ 10~ 29 gcm~ 3 , their deviation from 



spherical symmetry would result in a violation of the inequality (15). The situation 
worsens with superclusters of densities (10~ 20 — 10~ 31 )gcm~ 3 which are known to exist 

0- 

For clusters the answer to the question whether they are in gravitational equi- 
librium is less certain. With a typical cluster density of 10~ 18 gcm -3 , the crucial 
ingredient now is how much their shapes deviate from spherical symmetry. We spec- 
ulate that low density clusters with large off-sphericity could exist and a systematic 
survey (which goes beyond the scope of this letter) would be welcome. If most clus- 
ters have a geometry of low eccentricities and sufficiently large densities, the line of 
argument is similar to the case of elliptic galaxies discussed below. 

Elliptic galaxies with highest eccentricities are known to be of the type E7. 
This corresponds to (c/a) = 3.3. Why do we not see galaxies with a higher ratio than 
that corresponding to E7? Without invoking galactic dynamics, the answer could be 
the inequality (flBj) . In view of our results objects trying to reach the equilibrium, 
can do that, given enough time, by lowering the ratio of their axes. Apparently, 
galaxies have done exactly this. It should not be forgotten that our results include 
a positive cosmological constant A. If A = 0, from the point of view of the virial 
theorem, galaxies with a high ratio (c/a) could, in principle, be also in gravitational 
equilibrium. 

The significance of all results, presented in the paper, will grow with increasing 
value of A. In so-called quintessence models A is effectively a parameter which depends 
on the epoch ||. As such it could have been larger in the past. The investigation of 
structure formation in such models could be of relevance. 

Two issues remain to be explored in more detail in the future. Looking back 
at equation (|10D and assuming non-equilibrium, additional insight could be gained if 
we could solve the differential equation for Ijk(t). This requires that we succeed in 
obtaining more information on Kjk(t) and Wjk(t) in terms of new differential equa- 
tions or otherwise. Secondly, the expansion of the universe which can, in principle, 
oppose the gravitational steady state is in our equations only due to the A-term. It 
would be certainly worthwhile to make an attempt to include the full expansion. We 
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expect that this would strengthen our results. 
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